SPECTRUM OF THE VOLUME INTEGRAL OPERATOR OF 
ELECTROMAGNETIC SCATTERING 
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Abstract. Spectrum of the volume integral operator of the three-dimensional electromagnetic 
scattering is analyzed. The operator has both continuous essential spectrum, which dominates at 
lower frequencies, and discrete eigenvalues, which spread out at higher ones. The explicit expression 
of the operator's symbol provides exact outline of the essential spectrum for any inhomogeneous 
anisotropic scatterer with Holder continuous constitutive parameters. Geometrical bounds on the 
location of discrete eigenvalues are derived for various physical scenarios. Numerical experiments 
demonstrate good agreement between the predicted spectrum of the operator and the eigenvalues of 
its discretized version. 
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1. Introduction. Knowledge of the spectrum of an integral operator is impor- 
tant for the efficient numerical solution of integral equations. For example, it may 
help to choose an iterative algorithm, estimate its convergence rate, and construct 
a proper preconditioncr. The spectra of electromagnetic integral operators may also 
find direct application in physics and engineering: plasma studies, effective medium 
theories, and target identification. 

Numerical estimation of the operator's spectrum after discretization is not too 
practical as it requires even more computational resources than the complete numeri- 
cal solution of a scattering problem. Not to mention that such a spectrum will be that 
of a matrix, which only approximates the original operator. In these circumstances 
any useful analytical prediction is extremely valuable. 

In the pioneering papers £Q and [2] by R. E. Kleinman et. al. the scalar one-, 
two-, and (acoustical) three-dimensional cases, and the surface (boundary) integral 
equations were considered. Here we look at the integral equation arising in vectorial 
three-dimensional electromagnetic scattering on penetrable objects known as volume 
or domain integral equation. 

First results concerning the spectrum of the three-dimensional electromagnetic 
scattering problems were derived from the analytical Mie solution for a spherical 
scatterer. This case was analyzed in considerable detail in 0], and the so- 
called Mie-resonanccs were identified as eigenvalues. In 6 J. Rahola has compared the 
theory of the Mie resonances with the numerical spectrum of the volume integral oper- 
ator. Numerical experiments indicated, though, that apart from the expected discrete 
Mie resonances, the spectrum of the discretized volume integral operator contained 
eigenvalues clustered along distinct line segments in the complex plane. However, the 
ad hoc explanation of this phenomenon presented in |H] cannot be extended on the 
case of a general inhomogeneous scatterer. The latter requires a different and more 
rigorous approach. 

Our goal is to establish the connection between the operator's spectrum and the 
physical properties of a scatterer. Using the symbol of the electromagnetic volume 
integral operator we show that in addition to discrete eigenvalues the spectrum has an 
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essential continuous part as well. Fortunately, the symbol can be derived explicitly 
allowing for exact predictions of the essential spectrum for any Holder continuous 
anisotropic scatterer. In the isotropic case the essential spectrum is just the set of 
values admitted by the constitutive parameters of the object in R 3 . The approximate 
location of the discrete eigenvalues is estimated in the spirit of R. E. Kleinman et. 
al., and depends mostly on the extent and shape of the object. 

The paper is organizer as follows. In the next two sections we recall several more 
or less well-known facts from the theory of Fredholm operators with pointwise multi- 
pliers and simple singular integral operators. Driven by the electromagnetic problem 
at hand in Section 4 we consider 'combined' operators consisting of both pointwise 
multiplication operators and simple singular integral operators. Theory of such com- 
bined operators was developed by Mikhlin et. al. 0, jS] and we hope that our brief 
introduction will help to popularise their approach. The main purpose of these three 
sections is to introduce the concept of a symbol of a combined operator, which helps 
to establish the essential part of the operator's spectrum. Having mustered this tool 
we move on to the electromagnetic case where we construct the symbol of the volume 
integral operator, analyze the essential spectrum, and provide geometrical bounds on 
the distribution of discrete eigenvalues. Finally, we present some numerical exper- 
iments with the matrix version of the operator, and try to comprehend similarities 
and differences between the continuous and the discrete worlds. 

2. Spectrum of Fredholm operators with pointwise multipliers. The vol- 
ume scattering in one and two dimensions is described by integral equations with 
manifestly Fredholm operators, i.e. an identity operator plus a compact integral op- 
erator. This is also the case in the acoustic three-dimensional scattering and has to do 
with the fact that the singularity in the kernel of all these integral operators is weak, 
i.e. its order is smaller than three (dimension of the spatial manifold), whereas, the 
support of the scatterer is finite. For integral equations with Fredholm operators the 
questions of existence and uniqueness are linked, so that it is sufficient to establish 
the uniqueness of the solution only. 

With operators which are not of the form 'identity plus compact' the Fredholm 
property is established by proving the existence of a regularizer - an operator which, 
when applied to the operator in question, reduces it to the form 'identity plus com- 
pact'. There may be a left and a right regularizer, naturally. First, we clarify one 
important thing: the spectra of the original Fredholm operator and its regularized 
version are, in general, different. Indeed, consider the following Fredholm integral 
equation: 



where instead of the identity operator we have a (pointwise) multiplication with some 
continuous function a(x) ^ 0, x G D. The integral operator is presumed to be 
compact. In operator notation we shall write 




(2.1) 



x.x'GD 



[d(a) + K] u = v , 



(2.2) 



where d(a) indicates that we have a multiplication operator, analogous to a diago- 
nal matrix. The reqularized Fredholm equation is obtained from <|2.1[l by a trivial 
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operation of (pointwise) multiplication by a(x) , i.e. 

u(x) + a(x)- 1 J K(x,x')u(x') dx' = a(x)" 1 w(x). (2.3) 
x,x'eD 

Hence, the regularizer in this case is simply d(a) , as it gives the 'identity plus 
compact' form 

d(a)- 1 [d(o) + K] = I + diaj^K = 1 + K ', (2.4) 

where K' is a compact operator. The spectra of the two (original and regularized) 
operators are quite different though. Indeed, consider operator 

d(a)-XI + K, (2.5) 

and let A = a(xo), for some arbitrary xo € D. If A does not belong to the spectrum, 
then there exists a bounded inverse of (|2.5|l for that particular A. This means that if 
we consider equality 

[d(a) - XI + K] u = v , (2.6) 

then for all possible perturbations Av of the right-hand side and the corresponding 
perturbations Au of the solution there exists a constant C < oo such that 

||Au||<C||A«||. (2.7) 

Due to the contunuity of a(x) in D we can always choose a neighborhood of xo, say 
S(x ), such that \a(x) — a(x )| < e, for x 6 <5(x ). Consider perturbations Ait(x) 
which are zero outside 5(xq) and have a fixed norm ||Att|| = R, independent of S(xo). 
This can be a function of the following type: 

Att(x) ~ e 4X - k( H x e 6(x ). (2.8) 

or a similar oscillating function with varying spatial frequency k((5, e). Moreover, 
since K is compact, we may assume that our perturbations give \\K Au\\ < Ke, for all 
sufficiently small S(xo). Further, we have ||[a(x) — a(xo)]Au|| < Re. Hence, for the 
perturbation of the right-hand side we obtain 

|[A«|| < || (a - A)Am|| + \\KAu\\ <{R + n)e , (2.9) 

whereas, ||Au|| = R for all Au. For the condition (|2.7|l to hold we must be able to 
provide C < oo such that 

R<C(R+K)e, (2.10) 

for all e, and this is obviously not possible, since we can make e as small as we like 
by reducing <J(xo), while keeping R constant. Hence, the inverse of operator (|2.5|) is 
unbounded and A = ci(xq) belongs to the spectrum of operator [d(a) + K\. 

If continuous function a(x) is not constant on D, then the set of its values is an 
extended dense subset of the complex plane. Such a spectrum is called continuous 
as opposed to the usual discrete eigenvalues. In addition to continuous spectrum 
operator [d(a) + K] may have discrete eigenvalues as well. The regularized operator 
[/ + K'], however, is known to have discrete spectrum only. This follows from the 
fact that the spectrum of K' is purely discrete Thus, regularization destroys 
continuous spectrum. Finally, we notice that operator l|2.5|) cannot be regularized at 
the points of continuous spectrum. 
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3. Spectrum of simple singular integral operators. The operator of the 
three-dimensional electromagnetic scattering has the following form 



Au = d{a)u + A s d(b)u + Ku , (3.1) 

where K is compact and 

A s d(b)u = lim [ ffiWwx') dx' , (3.2) 

e-+0 J r 6 
x'GR 3 \<5(e) 



with r = x' — x|, and = (x'— x)/r. The principal- value integral above is understood 
as an integral over entire K 3 excluding a ball around x whose radius tends to zero in 
the limit. This integral exists under certain assumptions on 6(x)m(x) and /(©)■ The 
necessary condition is the following one: 



J /(0)dS = O, (3.3) 



where integration is performed over a unit sphere S. The Holder-continuity of &(x)it(x) 
in R 3 is another (sufficient) condition. Operator A s in l|3.2|l is called the simple 
singular integral operator. Its kernel is not an integrable function, and the best way 
to comprehend the meaning of such a p. v. -integral is to take the three-dimensional 
Fourier transform of A s u. Like this 



5(k) = lim / / ^5^ii(x')e- ?;k x dx'dx 
e-+0 J J r fi 

xel 3 x'eR 3 \<5(e) 

= lim J J ^ffl U (x')e- ik ( x - x ''e- lkx 'dx'dx 



xSK 3 x'GR 3 \<5(e) (^'^) 

lim J ^"^ e'^dy j u(x')e- ,;kx 'dx' 



ye 

= F(k)fi(k) . 



In other words, we have arrived at a familiar k-domain expression for the Fourier 
transform of a convolution. The function F(k) can only be found, if one knows the 
explicit form of /(©)• However, there is something to say even in the general case. 
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Consider this 



F(k) = lim / H-^l e -^y dy 
e^O J r 



y£K 3 \5( £ ) 

oo oo , 

lim J J ^ ~ e- irpcos7 dSdr = Urn / / /(-©) - * ^ dSdt 

e S pe S 

1 oo 

r r p —it cos 7 r r p—it cos 7 

limy j /(-©) dSd* + / / /(-©) dSd* 

pc s is 
00 1 (3-5) 

— it cos 7 p 1 

/(-0)^- dSd* + lim / / /(-©)- [ e - 4 * COS7 - 1] dS&t 

IS pe S 



lim J id* /" /(-e)dS 

pe s 



cos 7 C C ~\ 

f(-®)— t dSdt+ I J /(-©)- [e- 4 * co ^-l] dSd* 

IS OS 

Above we have itroduced spherical coordinate systems for both y and k, where p = k| , 
and used condition 1)3.3(1 . As one can see, since e — > tends to zero independently of 
p, i^(k) does not depend on p and, hence, it is independent of |k|. In other words, if 
F(k) depends on k at all, then it is only the direction of k, not its magnitude, i.e. 

F(k)=F(^j=F(&). (3.6) 

Note that this is only true if the order of singularity in 13.2|l is exactly three. If we 
had a weakly singular integral, then F(k) would depend on |k| explicitly. 

Let us now construct a regularizer for a simple singular operator ^4 S . Denoting 
by T and J 7 ^ 1 , correspondingly the forward and the inverse Fourier transforms, and 
by F(©) _1 the inverse of F(&), we can do, for example, the following: 

J r ~ 1 F(®)~ 1 J 7 A s u = ^ 1 F(©)- 1 J F 1 (0)J r u = T- x Tu = u . (3.7) 

In this way the original operator A s is reduced by the regularizer J 7 ~ 1 F(®)~ 1 T to 
an identity operator, i.e. the regularizer is also an inverse of the operator. 

To recover the spectrum of a simple singular operator we need to consider 



[A s -XI]u= lim [ ^Mu(x') dx' - Au(x) . 
e^O J r A 



(3.8) 



Repeating calculations l|3.4|) we arrive at the following equivalent version of the spec- 
tral problem: 

F(0)fi(k) = Au(k) , (3.9) 



G 



N. V. BUDKO AND A. B. SAMOKHIN 



Hence, generalizing this to matrix- valued F(&), we conclude that the domain s(A) 
in the complex plane consisting of all eigenvalues A of matrix F(&) for all on 
the unit sphere belongs to the spectrum of the simple singular integral operator A B . 
Note, although matrix eigenvalues are always discrete, they may vary continuously 
with 0, thus creating a continuous spectrum. Again, the operator of i|3.8|) cannot be 
regularized for A G s(A). 

4. Essential spectrum of combined operators. In the previous two sec- 
tions we saw that the spectrum and the regularization of operators are tightly linked. 
Namely, operators of the spectral problems (|2.5[) and l|3.8|l could not be regularized 
for A 6 s(A). This ia also true, if we consider the spectral problem for operator l|3.1|) . 
and, perhaps, for any operator. However, a rigorous proof of this link for a general 
case is too abstract, and we shall not provide it here. Taking it for granted we shall 
instead concentrate on the construction of a general regularizer. 

We have shown that the Fredholm operator with free term d(a) could be reg- 
ularized by d(a) _1 , whereas a simple singular integral operator is regularized by 
F~ 1 F(®)~ 1 T . However, neither of these approaches works with the general operator 
(|3.1I) arising in the electromagnetic scattering problem. Since this general operator 
is a combination of multiplication and simple singular operators (further - combined 
operator) we expect that a proper combination of the two approaches may work. 

Technique presented below is a simplified version of the approach suggested by 
S. G. Mikhlin [7], First thing to note is that, while searching for a regularizer, the 
actual form of the compact operator in the resulting 'identity plus compact' expression 
is not important. Secondly, a regularizer cannot be a compact operator, as the product 
of a compact operator with any other bounded linear operator is a compact operator. 
Finally, the product of any two combined operators is also a combined operator (plus 
a compact one, in general). For example the product of two operators of type (|3.1[1 
will look like 

A X A 2 = [d(oi) + A;d(6i) + Xi] [d(oa) + A\d(b 2 ) + K 2 ] 

= d(oi)d(a a ) + d( ai )A s 2 d(6 2 ) + A s 1 d(6i)d(a 2 ) + AJdC&O^dfa) + K , ( ' ' 

where all compact operators are hidden here in the generic compact operator K. 

To every combined operator there corresponds a function called symbol con- 
structed in accordance with the following set of rules. 

Definition 4.1. 

1. The symbol of a compact operator is zero. 

2. The symbol of a pointwise multiplication operator with a Holder continuous 
multiplier is the multiplier itself. 

3. The symbol of a simple singular integral operator is its Fourier transform. 

4- The symbol of the sum (product) of the aforementioned operators is the sum 
(product) of their symbols (the order of symbols in a product corresponding to 
the order of operators) . 
For the last property to hold it is sufficient to ensure that the multipliers are 
Holder-continuous functions in M. 3 , so that the products of simple singular integral 
operators and multiplication operators and the products of their symbols are well 
defined. According to Definition 14.11 the symbol of the combined operator (|3.1fl is 
just this function 



Symb [A] = $(x, 0) = o(x) + F(0)6(x) , 



(4.2) 
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and the symbol of the product l|4.1|l of two combined operators is 

Symb [A X A 2 ] = ai(x)a 2 (x) + a 1 (x)F 2 (0)6 2 (x) + F 1 (0)6 1 (x)a 2 (x) 

+ F 1 (0)6i(x)F 2 (0)6 2 (x) = $i(x, e)$ 2 (x, 0) (4.3) 
= Symb [Ai] Symb L4 2 ] . 

For our purposes the most important property of a symbol is that the correspon- 
dence between a combined operator and its symbol is one-to-one up to a compact 
operator. This means that to each combined operator there corresponds only one 
symbol, whereas to each symbol there may correspond more than one operator, but 
their difference would be a compact operator. 

Suppose now that we want to find a regularizer A r of a combined operator A, i.e. 

A r A = A 1 [d(a) + A s d(b) + K]=I + K' . (4.4) 

The symbol of this product can be constructed using Definition 14. II and is given by 

Symb [A 1 A] = Symb [A 1 ] |a(x) + F(0)6(x)l = 1 . (4.5) 

To this symbol there corresponds an operator of the form 'identity plus compact', i.e. 
should we succeed in obtaining such a unit symbol - the regularization is accomplished. 
Therefore, the symbol of a regularizer, say < i ) r(x, 0), must be an inverse of the symbol 
of the combined operator in question. Since the latter symbol may be a matrix- valued 
function, one may have to find an inverse of a matrix. 

In this simplified version of the Mikhlin's approach summarized in Definition 14.11 
we have limited ourselves to combined operators. Therefore, in order to be able to 
apply the multiplication property in 14.5|l we have to be sure that the regularizer, if 
it exists, belongs to the class of combined operators as well. This would be the case 
if the inverse of the symbol of a combined operator - [a(x) + F(0)6(x)] _1 - had the 
form of a symbol of a combined operator itself, i.e. contained only multipliers and 
Fourier transforms of simple singular integral operators. Fortunately, it follows from 
the Cayley-Hamilton theorem for matrices that the inverse of a nonsingular matrix can 
always be represented as a polynomial of this matrix. Obviously, any such polynomial 
contains only the products of the multipliers and the Fourier transforms mentioned 
above. 

Conditions on the existence of a regularizer thus become conditions on the ex- 
istence and boundedness of the inverse of the symbol of a combined operator for all 
x e K 3 and all on a unit sphere. 

The symbol will help us to recover the essential spectrum of combined operator 
A given by (2H) - Since operator A — XI is also a combined operator, its symbol is 
given by 

$(x, 0, A) = o(x) + F(0)6(x) - AI = $(x, 0) - AI , (4.6) 

where I is an identity matrix of the corresponding dimensions. This symbol is a 
parametric function of A as well as x and 0. Suppose that A coinsides with the 
eigenvalue of the symbol matrix ^(x, 0) for some x and 0. Then, operator A — XI 
cannot be regularized and, hence, does not have a bounded inverse. In summary 
(generalizing to matrix- valued symbols): domain s(A) in the complex plane consisting 
of all eigenvalues of symbol $(x, 0) for all x £ R 3 and all € § belongs to the 
(essential) spectrum of the corresponding combined operator A. If the eigenvalues of 
$(x, 0) vary continuously in the complex plane with x and varying correspondingly 
in M 3 and on S, then the essential spectrum is continuous. 
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5. The electromagnetic case. Now we turn to the electromagnetic three- 
dimensional case. In the frequency domain with time-factor exp(— iwt) the total 
electric field E in the presence of an inhomogeneous anisotropic electric-type scat- 
terer satisfies the following singular integral equation (see e.g. |12p 



E in (x,c) = 



I+~x(x,w) 



E(x,w) 

lim / G Q (x-x')x(x',w)E(x',w) dx' 

e^O J 

x'GK 3 \|x-x'|<e 

J Gi(x-x',w)x(x',w)E(x',o;) dx' 



where E m is a given incident electric field, and x is t ne tensor of relative electric 
contrast of the object with respect to a homogeneous isotropic background. Explicitly 
it is given by 

X (x,w) = — Lrrtx.w)-!, (5.2) 

where I denotes the 3x3 identity matrix and the constitutive parameters of the scat- 
terer are described by the tensor of the transverse medium admittance (per length), 
which in its turn contains the more familiar tensors of conductivity c(x, tu) and di- 
electric permittivity e(x, ui), namely, 

ry(x, bS) = <r(x, u) — iwe(x, ui) . (5-3) 

The isotropic background medium has parameters %(w) = <y\,{ui) — iue^uj) and 
Cb(w) = — iu) /j,b(u>) , where /Zb is the magnetic permeability. The so-called wavenumber 
kb(u>) of this medium can be found from k^(uj) = — T)(u>)£(uj) by taking a proper branch 
of a square root. Although the contrast function is different from zero only within 
a finite spatial domain D, we have extended the limits of integration in l|5.1|l to the 
whole of K 3 . This formal extention allows direct application of the results of the 
previous section. 

The kernels of the intergal operators in l|5.1|l contain tensors Go and Gi, the first 
of which is given by 

Mx) = 4^xP [3Q ~ I] ' (5 ' 4) 

where tensor Q = 0(0- ) can also be written in matrix notation as Q = 00 T . This 
tensor is an orthogonal projector, i.e. Q 2 = 00 T 00 T = Q. The second kernel 
contains tensor 

Gi (x, w) = [-a(x, uj) - 3/3(x, w) + 3 7 (x, u)} Q + [a(x, w) + /3(x, to) - 7 (x, u)] I, 

(5.5) 

with 

gife b (w)|x| e ife b (")|x| p ikb(u))\x.\ 



a(x,w) = k h (w) , | | , P(x,u) =ik h (cj) , 7(x,w) 



4tt|x| ' ^ v ' ' DV ' 4tt|x| 2 ' ]K 1 1 4tt|x| 3 ' 

(5.6) 
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In operator notation of the previous section our equation is written as 

= u in . (5.7) 



The operator of this equation acts in the Hilbert space L^(M 3 ) of square-integrable 
vector- valued functions. This choice of the functional space is physically motivated 
by the energy considerations. All components of the matrix-valued function x(x, u>) 
are presumed to be Holder-continuous in R 3 , i.e. not only inside the object, but also 
across its outer boundary. As previously K denotes a compact operator, which in this 
case is the last integral operator in (|5.1I) with weakly singular tensors H5.5[) - I|5.6[l in 
its kernel. Although, in general, weakly singular integral operators are not compact 
on K 3 , the kernel of the last operator in (|5.1(l contains the contrast function as well, 
which effectively restricts the spatial support of the domain of the integral operator 
to D. Therefore, the total operator is compact. 

There exist other formulations of this problem. For example, for isotropic scat- 
terers with differentiable constitutive parameters the regularization of the integral 
equation can be carried out explicitly, see HH|, However, from the discussion 

of the previous sections we know that the spectra of the original and the regularized 
operator may be very different. 

To find the spectrum of our operator we first need to derive its symbol. The most 
laborous task here is, obviously, the Fourier transform of the simple singular integral 
operator. The function /(©) required in (|3.5|) is in our case 

/(©) = i-[3Q-I] . (5.8) 

It is a matrix-valued function and therefore its Fourier image F(®) is also a matrix- 
valued function. The complete symbol of the electromagnetic volume integral operator 
was first obtained in using the technique suggested in [5], where /(©) is expanded 
in terms of spherical harmonics. Comparing the result of |12| with the present formu- 
lation we find that 

F(0) = ~I-Q, (5.9) 

where Q = 0(0- ), or in matrix notation Q = 00 T . Now using Definition 14. II and 
expression (|5.2() we obtain the symbol of the electromagnetic volume integral operator 
for an anisotropic scatterer as 

*(x,e) =I + Qx(x,o;) =I-Q + Qffr(x,w), (5.10) 

where both \ arL d fjr — T)/Vb are matrix- valued functions and in general do not com- 
mute with matrix Q. The inverse of this symbol can be found explicitly as 

$- 1 (x,0)=I + r 1 Q-r 1 Q'?r(x,a;) ) (5.11) 

where 

e(x,0) = T 7 7l .(x, W )0. (5.12) 
In isotropic case, where r/(x, lo) is a scalar, the inverse of the symbol looks even simpler: 
^otro Pi c( x >®) = I - Q + Qv7 1 ( x > iJ ) = I + Qx'(^oj) . (5.13) 
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with x' given by 



X'(x,0,) = ^\-1. (5.14) 
?7(x, u>) 



Although it is not our primary concern, now we can find the explicit electromagnetic 
regularizer. Notice that expression Ij5.13|l is identical in form to the symbol Ij5.1t)|l 
of the original operator. This means that the operator corresponding to $ _1 has 
the same form as the operator corresponding to $. Hence in the isotropic case the 
regularizer is 



isotropic = d(I + -x') + A s d( X ') , (5.15) 



with x' given by (|5.14|) . 

The essential spectrum consists of the eigenvalues of the symbol matrix $(x, 0) 
for all x £ M 3 and all on the unit sphere §. Direct computaions show that in the 
anisotropic case these eigenvalues are Ai = A2 = 1, and A3 = T ^ r (x, w)0, and they 
depend on both x and 0. Whereas in the isotropic case, where r] T is a scalar, the 
eigenvalues are Ai = A 2 = 1 and A3 = ?7 r (x, lv), and they depend on x only. Hence in 
general the essential spectrum is given by: 

Af 8 = X C 2 SS = 1 , Xf s = © T r/ r (x, oj)& , x e R 3 , e § . (5.16) 

Notice that there is also a parameteric dependence on the angular frequency u> in both 
isotropic and anisotropic cases, which matters if the scatterer is conductive and/or 
dispersive. 

The object parameters are presumed to be Holder continuous functions of x on 
R 3 . This means that even with a homogeneous scatterer, characterized by 7y r (x, w) = 
const, x £ D, there exists a smooth transition between r\ x and 1 the background 
medium parameter. Then, the essential spectrum must connect 1 and point r/ r in the 
complex plane. Morover, for any inhomogeneous scatterer the essential spectrum will 
emerge from point 1 on the real axis, and continuously span all values admitted by 
T 7? r (x,w)0 for x 6 D and 8e§. 

Apart from the essential spectrum our operator may have discrete eigenvalues, 
which correspond to nontrivial solutions of the problem 



d(l+^x)+A s d( X )+K 



u-Xu = 0. (5.17) 



We shall consider here only the isotropic case. In this respect it is convenient to use 
the standard sufficient uniqueness condition for the original problem (|5.7|l . which in 
terms of the admittance function ?y(x, to) given by (|5.3|) can be written as 

Re[7j(x,o;)] > 0, xefl. (5.18) 

Under this condition the null-space of the volume integral operator (operator in square 
brackets in (|5. 17(1 ^1 is known to be trivial. Considering the case where A ^ 1 we can 
re-write Ij5.17|l as 



d(I+±x')+A s d(x')+K> 



u = , (5.19) 
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where the only difference with the original operator is the new contrast function x' 
given by 

x(x,w,A) = - — = 1 = — 1. (5.20) 

1 — A 7/b(l-A) Vb(U) 

And since the eigenvalue problem has been reduced to the problem about the non- 
trivial null-space of the original operator we can use an 'inverse' of the condition 
(|5~TS|> . i.e. 

Re [r/(x, oj, A)] < 0, xefl. (5.21) 
From here we derive the following bound on the eigenvalues: 
Re7y(x, w) — [Re?7(x, u) + Rer?b(w)] ReA 

„ (5.22) 
- [Imr?(x,w) -Imr? b (w)]IniA + Rer;b(w)|A| 2 < 0, xefl. 

In addition we note that 

|A| < \\A\\ < oo, (5.23) 

since A is bounded, and that A = is not an eigenvalue if the scattering problem has 
a unique solution (which we presume from now on). 

The only natural infinite homogeneous background medium is vacuum, where 
(Tb(w) = and £b(w) = £o- hi particular cases, such as scattering from subsurface 
objects, other homogeneous backgrounds can be used as an approximation. We shall 
deal with those later. Let us first consider the natural case, and in addition presume 
that cr(x, lu) = er(x) and e(x, lu) ~ e(x), and are real-valued functions. Further we 
also assume that the contrast in permittivity e(x) — £b is positive for all x 6 D. These 
assumptions hold for a static model of a conducting dielectric, which is applicable at 
low frequencies where the dispersion of constitutive paremeters is insignificant. In 
this case the eigenvalue bound (|5.22|) reduces to 

<j(x) - o-(x)RcA + ui [e(x) - e ] ImA < 0, xefl. (5.24) 

In Fig. 5.1 and Fig. 5.2 (left) we present geometric estimates derived from this formula 
for the following three situations: 

e (V) 

a(x) =0, 1< — < M ; (5.25) 

cr(x) = const > 0, e(x) = const > e ; (5.26) 
cr(x) > 0, e(x) > e ; (5.27) 

where xefl. The other parameters in these figures are given by 

M= — +i—, N = —i — - — — -, (5.28) 

£o ^£o cj(e-eo) 

c(x) 

Ni = -i max = , \ ' -, (5.29) 

xec oj [e(x) - e ] 

N 2 = -i min -, (5.30) 
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Fig. 5.1. Location of the spectrum for a homogeneous object in vacuum. Left: nonconducting 
dielectric. Right: conducting dielectric. 



Fat solid lines and curves schematically outline the essential part of the spectrum, 
which in the case of a homogeneous object we presume to be straight lines connecting 
the real 1 and the corresponding value of rj/rfb- As one can see, in the natural (vacuum) 
background medium, the bounds are wedge shaped, same as in pp. In the case of an 
object with losses these bounds safely separate the domain of eigenvalues from the 
zero of the complex plane. This is important since the uniqueness condition can only 
guarantee that there are no eigenvalues equal to zero, but does not tell if there are 
any eigenvalues close to zero, while such eigenvalues may cause instabilities in the 
numerical solution and slow down the convegrence of an iterative algorithm. From 
this point of view the worst case is the lossless object or an inhomogeneous object 
with lossless parts. 

Now let us analyze the popular approximation, where the homogeneous back- 
ground medium is considered to be conducting, i.e. <7b > and/or has arbitrary 
constant permittivity £b > £o- Such approximations are often used if a scatterer is 
located inside a large but finite and more or less homogeneous host medium, say, 
inside the Earth or a human body, and the reflections from the outer boundary of the 
host medium can be neglected. We again consider the static model of a dielectric and 
from (|5.22(1 we derive the following estimate 



(5.31) 



which shows that all eigenvalues are situated inside a circle. Parameters of this circle's 
center are 





2 




2 




Re A — ^ 


+ 


ImA+f 


< 


(f) 













A =^ + 1, B = — [e(x) - e b ] , 

Ob fb 



and its diameter can be found from 

>(x) 



D z 



1 2 



o-b 



(5.32) 



(5.33) 
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Figure 5.2 (right) corresponds to a particular situation where a nonconducting object 
(cr(x) = 0) with no contrast in permittivity (e(x) = Sh) is immersed in a conducting 
background. This can represent an air pocket or a low-contrast plastic landmine in 
wet soil. The constant M2 is given by 

M 2 = JL = ^!fL_ _ (5.34) 

On one hand, the spectrum is now explicitly bounded allowing, in particular, for a 
more precise convergence estimate when an iterative method is employed to solve 
the scattering problem. On the other hand, the circle extends towards the zero of 
the complex plane as in Fig. 5.2 (right), meaning that, although there is no zero 
eigenvalue in the spectrum, some of the eigenvalues may get close to zero, and render 
the scattering problem unstable. 

The last point to note is about the relative weight of the compact and singular 
parts in the operator of 1)5. lfl as a function of angular frequency oj. The coefficients 
a and (3 in (|5.ti|l depend on the wavenumber fcb(w), which in its turn is proportional 
to oj. Hence, the kernel of the compact operator gets more 'weight' as frequency 
increases, and the norm of the operator increases as well. As the compact operator 
delivers only eigenvalues, we can expect that the eigenvalues will spread out at higher 
frequencies. 

6. Numerical experiments. There is only a limited correspondence between 
the spectra of an integral operator and of its discretized (matrix) version. First of 
all, matrices do not have continuous spectra, but only discrete eigenvalues. Hence, we 
should not expect to see the lines and curves of the previous section. On the other 
hand, the continuous spectrum may serve as an accumulation area for the eigenvalues, 
as the latter should converge to the spectrum of an operator in the continuous limit. 
Let us first see if it is indeed so. 

We discretize the operator of (|5.1|l using the standard collocation method (see 
e.g. [B], |12|). which gives an order h 2 accuracy, where h is the size of an elementary 
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Fig. 6.1. Different discretization (and quantization) levels for three cubes (homogeneous, two-, 
and three- layered). The side length of all cubes is a f» A mo( j/20. Inhomogeneity becomes 'visible' 
only at an appropriate discretization level. 

cubic cell. The scatterer is discretized on an TV-point homogeneous grid, and the 
resulting matrix is 37V x 37V and is completely filled with complex numbers. The 
matrix eigenvalues are computed using the standard Matlab function eig. 

According to our theoretical predictions the essential spectrum, continuous or 
not, dominates at lower frequencies (small objects). Therefore, we shall first consider 
an object whose extent is sufficiently small. The background medium is vacuum 
and the frequency is set to u)/(2w) = 1 GHz. The scatterer has the shape of a 
cube with side length a w (A moc i/20), where A me d is the medium wavelength, see 
Fig. 6.1. According to the common ad hoc discretization rule such a small object 
can be modelled as a single cell (Fig. 6.1, left). This, of course, would only give an 
order h 2 approximation of the field at the geometrical center of the scatterer. In 
an attempt to compute the field at other points inside the object one may wish to 
introduce a finer discretization. Then, if the object is homogeneous, i.e. r/(x) = const, 
the value of the consitutive parameter would be the same at all internal grid points 
(Fig. 6.1, upper row). If, however, the object is inhomogeneous, then the matrix 
will include the new, refined, grid values of f?(x„), n = 1,...,N (Fig. 6.1, middle 
and lower rows). Obviously, there is a link between the spatial discretization and 
the quantization of constitutive parameters, and we have observed an interesting 
phenomenon related to these two processes in our numerical experiments. Namely, 
if we refine the spatial discretization while keeping constitutive parameters constant, 
then the (low-frequency) spectrum converges to a set of very dense line segments 
connecting the real unit and the corresponding values of rj/rjt,. In this way we seem to 
model an object consisting of one or more homogeneous parts without really telling 
what is the behavior of the constitutive parameter across the interfaces. If, on the 
other hand, the quantization of the constitutive parameter results in new values at a 
finer level of discretization, then the line segments in the spectrum are less pronounced. 
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Fig. 6.2. Effect of discretization (and quantization) refinement on the matrix eigenvalues in the 
low-frequency (small object) regime. Plots correspond to the objects shown in Fig. 6.1, constitutive 
parameters n/rj^ are given as circles. 



Consider, for example the case of three different cubes having the same size, and 
the following sequence of grid steps: h = a, a/2, a/4, a/8. The matrix dimensions are 
then, correspondingly: 3 x 3, 24 x 24, 192 x 192, and 384 x 384. The objects and the 
eigenvalues are presented, respectively, in Fig. 6.1 and Fig. 6.2. The upper row in both 
figures corresponds to the model of a completely homogeneous cube, the middle row 
models a cube consisting of two layers, and the bottom row corresponds to a three- 
layered cube. The layered structure becomes visible at h = a/2 (second column), 
and the difference between the two and the three layers is only seen at h = a/4 
(third column). At h = a/8 (last column) we had no more changes in the constitutive 
parameters of all three objects. The constitutive parameters rj/r]h are shown as circles 
in Fig. 6.2. The dense lines appear as soon as we stop refining the quantization, but 
keep refining the spatial discretization. Figure 6.3 illustrates the marginal influence 
of the object's geometry on the distribution of eigenvalues at low frequencies (same 
frequency as before). Both objects (top and bottom row) consist of two homogenous 
parts with constitutive parameters given by circles and have the same outer shape 
(long parallelepiped, roughly A/5 x A/42 x A/42). The difference between the objects 
is in the geometry of the parts. The upper row in Fig. 6.3 corresponds to the object 
consisting of two equal parts (roughly A/5 x A/42 x A/85 each), which divide the 
parallelepiped along its longer dimension. Whereas, the lower row corresponds to 
the object consisting of two unequal parts (roughly A/21 x A/42 x A/42 and A/7 x 
A/42 x A/42) dividing the parallelepiped across its longer dimension. From left to 
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Fig. 6.3. Effect of the object's shape on the matrix eigenvalues in the low-frequency (small 
object) regime: same constitutive parameters rj/r]^ (circles), but different shapes. From left to right: 
pictures of the objects, and eigenvalues for two discretization levels: h = A me< j/70 and A me< j/140. 
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Fig. 6.4. Effect of discretization on the matrix eigenvalues at higher frequencies. Object 
homogeneous cube. Left: discretization level h = A met j/5. Right: discretization level h = A me< j/10. 



right two discretization levels are presented: h = A mc( j/70 and h = A mc d/140. With 
both objects we notice the appearence of the line segments connecting the real unit 
and the values of r)/rfh. The difference seems to be among the other eigenvalues 
situated between the line segments. These and many other low-frequency numerical 
experiments suggest that the line segments observed here and in [S] are the matrix 
analogue of the operator's essential spectrum, and that they depend only on the grid 
values of r](x)/rib- 

Figure 6.4 gives the spectrum for a homogeneous cube at higher frequencies. The 
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Fig. 6.5. Spectrum in exotic cases. Left: vacuum cube in a lossy background medium. Right: 
cube with negative permittivity and small losses in vacuum background. 



r]/i% of the object is shown as a circle. The length of the cube's side is chosen 
to coinside with the medium wavelength in this case, i.e. a = A mo d- One can see 
that now we have not only the dense line segment as in Fig. 6.2 (top-right), but 
also a few off-line eigenvalues, which appear within the predicted bounds of Fig. 5.1 
(right). Two discretizations are shown: h — A mo d/5 (Fig. 6.4 left), and h = A me d/10 
(Fig. 6.4 right). Normally, neither of these would be considered a "proper" level of 
discretization. However, here we observe an interesting phenomenon related to the 
essential spectrum, which partly justifies a much more relaxed attitude of practitioners 
to the discretization of the electromagnetic volume integral equation, as opposed to 
the tough requirements on the discretization of differential Maxwell's equations via 
the finite-difference approach. With h = A me d/5 we have in total 375 eigenvalues 
most of which are clustered along the line segment of the essential spectrum. With 
h = A m ed/10 we already have 3000 eigenvalues, however, all 'new' 2625 eigenvalues 
keep filling the same line segment. Of course, the off-line eigenvalues do shift a little, 
but their number seems to remain the same for both discretization levels. This and 
similar experiments at higher (resonance) frequencies show that the refinement of 
discretization has no significant effect on the spectral radius of the matrix, which is 
what is expected from an integral equation formulation. 

We conclude our numerical experiments with the result corresponding to the 
conductive background case described in Fig. 5.2 (right). As one can see from Fig. 6.5 
(left) the matrix spectrum is, indeed, quite accurately described by the predicted 
circle. 

7. Conclusions and possible applications. The spectrum of the volume inte- 
gral operator of three-dimensional electromagnetic scattering contains both the essen- 
tial continuous part and discrete eigenvalues. The apparent difference with the one- 
and two-dimensional electromagnetics, where the spectrum is purely discrete, stems 
from the presence of a strong singularity in the kernel of the integral operator in the 
three-dimensional case. We have shown that the essential spectrum is explicitly given 
by (|5.16(l for any anizotropic scatterer with Holder continuous constitutive parame- 
ters. Knowing the spectrum one can find, for instance, the relaxation parameter which 
provides the optimal convergence of the over-relaxation iterative algorithm. Due to 
the fact that the well-described essential spectrum dominates at lower frequencies we 
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suggest that the over-relaxation method (which is very cheap from the computational 
viewpoint) should be used with quasi-static problems. 

However, there exists a much more exciting direct application of this knowledge 
as well. For some reason various electromagnetically 'exotic' media are at the core 
of the present day research. Take for example left-handed materials also known as 
media with negative refractive index. These materials are supposed to be highly 
dispersive, so that for a certain range of w's both the dielectric permittivity and the 
magcntic permeability happen to have negative real parts. Although we have not 
discussed the magnetic case here, our preliminary calculations indicate that all the 
present conclusions hold, and that the magnetic properties of the object result in 
another similar contribution to the essential spectrum. Then, the spectrum of the 
'left-handed' object will not only contain points in the left part of the complex plane, 
but also lines or curves connecting these points with the real unit. For a low loss 
material, which for obvious reasons is the ultimate goal of experimentalists, the lines 
of essential spectrum may then proceed dengerously close to the zero of the complex 
plane as shown in Fig. 6.5 (right), where e/eb = —1 and a = 0.001 S/m. From the 
physical viewpoint this would mean that the electromagnetic field is unstable in such 
media. Without losses {a = 0) the line would go right through the zero, meaning that 
the solution (electromagnetic field) does not exist at all. May be that is the reason 
why we do not observe many left-handed substances in Nature? Another interesting 
application where the explicit knowledge of the essential spectrum may be of help is 
magnetically confined fusion plasma. 
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